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Abstract 

The double series approximation method of Bonnor is a means for examining 
the gravitational radiation from an axisymmetric isolated source that undergoes 
a finite period of oscillation. It involves an expansion of the metric as a double 
Taylor series. Here we examine the integration procedure that is used to form 
an algorithmic solution to the field equations and point out the possibility of 
the expansion method breaking down and predicting a singularity along the axis 
of symmetry. We derive a condition on the solutions obtained by the double 
series method that must be satisfied to avoid this singularity. We then consider 
a source with only a quadrupole moment and verify that to fourth order in each 
of the expansion parameters, this condition is satisfied. This is a reassuring 
test of the consistency of the expansion procedure. We do, however, find that 
the imposition of this condition makes a physical interpretation of any but the 
lowest order solutions very difficult. The most obvious decomposition of the 
solution into a series of independent physical effects is shown not to be valid. 

1 Introduction. 

In 1959 Bonnor introduced the doubles series method as a means for looking at 
the gravitational radiation from an isolated source that oscillates for a finite period 
[0. This involved expanding the metric as a double power series in two parameters 
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and solving the vacuum Einstein field equations by successive approximations. This 
first approach dealt with isolated sources emitting spherical gravitational waves and 
showed that the source lost mass at a rate equal to that of energy radiated. Following 
the publication of the Bondi metric in 1960 0, a metric well suited to this problem, 
Bonnor and Rotenberg refined and generalised the double series method to the ax- 
isymmetric case 0. They were then able to reproduce the earlier result concerning 
mass loss and to show that the momentum generated in the source is equal and op- 
posite to that removed by the gravitational waves. They were also able to assert the 
existence of wave-tails in Bondi coordinates. These represent the backscattering of 
the gravitational radiation by the spacetime curvature induced by the source. The 
tails were interpreted as incoming radiation in a further paper on the double series 
method by Hunter and Rotenberg Q. Elsewhere, [||, I considered the interaction of 
these incoming wave-tails with the Schwarzschild source. I found that the wave-tails 
have no permanent effects on the source - all changes in the metric die off at least 
like u~ 2 (u being retarded time). I also showed that the incoming radiation does not 
alter the mass of the source itself. 

The double scries method itself has influenced the formalism of Blanchct, Damour 
and Iyer Q, in which the spacetime is divided into a near-zone region, a far-zone 
region and a region of overlap. In the far-zone a post-Minkowskian expansion is 
undertaken that is a generalisation of the double series method to more physical, 
non-axisymmetric sources. A review of the BDI formalism has recently been given by 
Blanchet @. 

In this paper I derive a consistency condition for the double series method. If this 
condition is not satisfied at a given level of approximation, a singularity appears in 
the solution along the axis of symmetry. While it might be possible to interpret this 
singularity as an infinite pipe or strut, this would go against our original criterion 
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imposed here, namely that the system be isolated. The occurrence of this singularity 
would therefore represent a breakdown of the double series method applied to isolated 
sources. The importance of the condition lies in the fact that at a given order of 
approximation its satisfaction depends not on the solution of the field equations to 
that order, but on the previously obtained lower order solutions. The condition is 
more a test of the double series method itself than of any particular solution obtained 
by the method. 

We then extend the work on the double series method using a power series package 
developed for SHEEP to carry out the otherwise prohibitive algebraic calculations . 
We verify that the consistency condition is satisfied at the third order and obtain the 
corresponding solution. We are then able to show that to fourth order the consistency 
condition is also satisfied. However, when considering specifically the wave-tail - 
wave-tail interaction we find that the condition fails. This component of the solution 
generates terms in the metric that are singular along the axis of symmetry. We are 
able to show that there is an exact cancellation that removes these singular terms in 
the full solution. This suggests that the double series method is not decomposable in 
the sense that it cannot be easily separated into a series of individual physical effects 
but that there is a complicated mixing of these effects occurring. This in turn leads 
to difficulties in the physical interpretation of the solutions beyond the lowest orders. 

The order of this paper is as follows. In section 2 we describe the source and the 
expansion of the metric coefficients that is characteristic of the double series method. 
In section 3 we outline the algorithm that is used to solve the linearised field equations 
and in section 4 we describe how this algorithm leads to the possible occurrence 
of singularities along the axis of symmetry. We then derive the condition for the 
avoidance of these singularities. Section 5 reviews previous work on the double series 
method and verifies that the solutions obtained satisfy this consistency condition. In 
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sections 6 and 7 we extend the solution, for a source with only a quadrupole moment, 
to the third and fourth orders and verify that these solutions are consistent with our 
assumption of an isolated source to fourth order. 

2 The Double Series Method. 

We consider an isolated axisymmetric source, such as two particles oscillating along an 
axis. The system is at rest for u < u±, vibrates smoothly for the period u\ < u < U2 
and is again at rest for u > ii2, though not necessarily in the same state as before u\. 
We are interested in the behaviour of the system at large radial distances and so do 
not specify the source in detail. We do not describe the mechanism that drives the 
oscillations since this should not have an effect on the asymptotic behaviour of the 
solution. 

It is assumed that there are two characteristic measurements of the system, m, 
the source mass, and a, a characteristic length (the separation of the two particles 
in the example mentioned above). The fundamental assumption of the double series 
method is that there exist solutions of the vacuum field equations that admit Taylor 
expansions in m and a about m = 0, a — 0. We do not offer a proof that this is the 
case, but the linear approximation suggests that this is a valid assumption to make. 

We assume that the system can be modelled by Bondi's radiative metric f| which 
can be written in the form: 

ds 2 = -r 2 (Bd9 2 + Csin 2 9d<p 2 ) + Ddu 2 + IFdrdu + 2rGd9du (2.1) 

with C = B^ 1 and where the metric coefficients are functions of u, r, 9 . Here we are 
using Bondi coordinates labelled (x\, X2, £3, £4) = (r, 9, <fi, u). 
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We make two further assumptions in order to start the approximation procedure. 
The first is that if we remove the mass (to = 0) our metric should reduce to flat 
spacetime in Bondi's coordinates: 

ds 2 = -r 2 {d8 2 + sin 2 9d(f> 2 ) + du 2 + 2drdu. (2.2) 

The second assumption is that if we set a = then the spacetime should reduce to 
Schwarzschild space in Bondi coordinates: 

2m 

ds 2 = -r 2 (d9 2 + sin 2 6d<j) 2 ) + (1 )du 2 + 2drdu. (2.3) 

r 

Our fundamental assumption leads to the following representation: 
B = L+XE 

p=l S=0 

c = i+xx c ™ p « s 

p= i s =o 

OO OO / . 



D = 1 + ^^ m p a s 



OO OO / ^ 

(ps) 



F = 1 + ^^ F m p a s 

p=l s=0 



OO OO , , 

(ps) 



G - XX G m p a s . (2.4) 

p=l s=0 

(10) 

To zeroth order in m this corresponds to the metric (|2.2| ). The choice D = — 2r~ 1 1 

(s0) (sO) (sO) (sO) (sO) 

with all other B , C , D , _F , G , (s > 0) vanishing, ensures that the metric reduces 



to (2^) to zeroth order in a. Here again the coefficients in the expansion (2A) are 



n • • (P s ) 

functions of u, r, 9. We now introduce the notation g a b to represent the terms in the 
metric proportional to m p a s . By the (ps) case we shall mean the solution of the field 
equations to order m p a s . Where there is ambiguity as to which products give rise to 
m p a s terms (for example in the (24) case we have 

ma 2 x ma 2 ,m x ma 4 , ma x ma 3 (2-5) 



terms), we shall refer to these products explicitly, e.g. the (12) x (12) case, the 
(10) x (14) case or the (11) x (13) case in the example above. 



Having made the expansion (2.4) above, the vacuum field equations are of the 
form: 

nth - (2.6) 

where the left hand side is linear in <f a b (the unknowns) and the right hand side is non- 

(qr) 

linear in o {, , q < p, r < s which are known from previous steps in the approximation. 
The explicit form of these equations is given in the appendix. 

We have now arrived at a system of linear partial differential equations. In the 
next section we describe the algorithm that is used to obtain solutions to this system. 

3 The Algorithm for Solving the Field Equations. 



It can be shown that the field equations (2.6), three of which are identically satisfied 



due to our imposition of axial symmetry, admit an algorithmic solution. The deriva- 
tion of this algorithm is given in the appendix. Here we give the bare details that are 
required for the derivation of the consistency condition. In order to ease the notation 
we consider the general m p a s case and omit the superscripts (ps). The algorithm 
consists of three steps: 
Step 1. 

Obtain F, the coefficient of gi2 of order m p a s , as defined by ( [2.4| ) by the relation- 
ship: 



F = -\J rHdr + r)(9,u) (3.1) 
where rj is a function of integration and H is one of the non-linear terms, i.e. one 



component of the right-hand side of equations (2.6). 
Step 2. 
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Solve for D the 'pseudo-wave' equation 



UD = D 11 -2D 14 + 2r- 1 (D 1 +D 4 )+r- 2 (D 2 2 + D 2 cot8) 
= -K + 2(F U + 2r- 1 F 4 ) + 



2r~ 



J{r 2 (N - 2F M ) + (F 22 + F 2 cot 9)}dr + x (t 



(3.2) 



where if, N are non-linear terms, x is a function of integration and F is defined by 



(3.1). The pseudo-wave operator denoted □ is similar to the usual d'Alembertian 
wave operator of flat space in the Bondi coordinates used here, the only difference 
coming in the sign of the 2D^r~ x term. 
Step 3. 

Obtain the remaining metric coefficients from the expressions 



G = r 1 cosec 



9 J sin9 J r 2 {N- 2F 1A )dr + r 2 D 1 + X 
/ F 2 dr + v(r, u) cosec# 



d6 



(3.3) 



B = cosec 2 9 J sin 2 - J {rL + 2r~ 1 (F 2 - G)}dr + F 2 -G-rG x 
+r(r, u) cosec 2 + n(0, u) 



d0 



(3.4) 



where is, r, \i are further functions of integration and L is another non-linear term. 

The five functions of integration that arise in the double series method, namely 
77,X,^,t and /z, are chosen to satisfy two criteria. Firstly they should ensure that 
the solution satisfies the full set of field equations (^1]) (the derivation of the above 
algorithm depends on only 4 of the 7 non-trivial field equations). Secondly they are 
used to impose certain boundary conditions on the system, namely that the metric 
be Minkowskian at spatial infinity, and that the metric remain non-singular on the 
axis of symmetry except at the origin. A sufficient condition for the latter is that 
B cosec 2 9, C cosec 2 9, D, F,G cosec 9 are of class C 2 near sm9 = 0. In the next 
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section we show that unless a particular condition is satisfied this criterion cannot be 
achieved. This imposes a consistency condition on the double series method. 



4 A Consistency Condition. 



We now examine the possibility of a singularity arising along the axis of symmetry 
during the integration of the field equations contained in the algorithm of the previous 
section. This singularity cannot be reconciled with our assumption of an isolated 
source and its existence would imply a lack of consistency in this application of the 
double series method. In later sections we will show that this condition is indeed 
satisfied to order m 4 a 4 . Here we derive the condition itself. As in the previous 
section we consider the general m p a s case and drop the superscripts (ps) to ease the 
notation. 



We suppose that we have obtained a solution F of (3.1) and a solution D of the 



psuedo-wave equation (3.2). Since the wave operator □ defined by (3.2) preserves 
the Legendre polynomials in cos#, it is natural to make the following expansion of 
the theta dependence of the solution: 

oo oo 

D = iD(u,r)Pi(cos6), F = ^ iF(u, r)Pj(cos0) 

t=0 i=0 
oo oo 

N = J2 iN(u,r)Pi(co S 6), X = P * ( cos 

i=0 i=0 

where Pi is the ith Legendre polynomial. 

We now proceed to step 3 of the algorithm and obtain an expression for G from 
(0): 

00 r r r 

G = ^V -1 cosec0 / sin0Pi(cos0) / r 2 (,N - 2 t F 14 )dr + r 2 l D 1 + lX d6 
i=o J U 

oo „ 

+ ^V- 1 P J (cos0)' / iFdr + v(u,r)cosec6, (4.2) 
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where here denotes differentiation with respect to thcta. Let us now take the first 
term [i = 0) of the series in the above expansion. This is 



r cot ( 



r 2 ( N -2 F 14 )dr + r 2 D 1 + aX 



which is singular at sinf? = 0, i.e. along the axis of symmetry. Since 



Pi (cos 9) sin 



(4.3) 



(4.4) 



contains a factor sm9 for i > 1, the expression (|4.3|) is the only singular term in the 



series contained in (4.2). Moreover, this singularity cannot be removed by a specific 



choice of the function of integration v in (4.2). The choice 



v = ±r 



r 2 ( N -2 F u )dr + r 2 D 1 + oX 



(4.5) 



removes the singularity along either the semi-axis 9 — 0, or the semi-axis 9 = n but 
not both. Thus we find that the double series method predicts a singularity along the 



axis of symmetry unless the term (4.3) vanishes, i.e. unless 

Ui = -r- 2 Jr 2 ( a N~2 Q F u )dr - r~ 2 oX . (4.6) 

We now take this expression for o-Di and substitute it into the pseudo-wave equation 
([T^) which leads us to an expression for qD^: 

D 4 = r[\ Q N -\ Q K] +2 F 4 . (4.7) 

Now using the fact that partial derivatives commute we can obtain two expressions 



for 0-D14 from (4.6) and (4.7) which must be equivalent. This leads us to the following 
condition that must be satisfied if a singularity along the axis is to be avoided: 

X ' = r2 [ I O-Wll ~ 5 oK\i + I 0-^44 — 5 0-^14 ] + 

r[2 Q N 1 -2 Q K 1 + oNi-loHi] +[ N- K] 
= 0. (4.8) 



This is a necessary condition for the avoidance of singularities at the (ps) level. It is 
also a sufficient condition for the avoidance of singularities in Q in the sense that 



if (4.S) is satisfied then a particular choice of the functions of integration x, V can 



(ps) 

always be made to ensure that there is no singularity along the axis in G ■ There is 
the additional possibility of the integration procedure leading to a singularity along 

(ps) 

the axis in B ■ The avoidance of this singularity leads to a condition on the P\ (cos 6) 
coefficients of the non-linear terms and on qD\. However, this condition turns out 



not to be as simple as (4.8), nor is it relevant to the rest of the work in this paper as 



we shall see in the next section. 

It can now be seen that the condition fll-SP depends only on the non-linear terms 

(qr) 

qN, qH, qK, i.e. only on the metric coefficients g a b, Q < P, t < s, that have 
been found by previous approximations. The occurrence of a singularity along the 
axis would, for the type of sources we have in mind here, appear unphysical. Any 
interpretation of this singularity as a strut or pipe extending along an infinite semi-axis 
would contradict our assumption the the source be isolated. Therefore in solving the 
(qr), q < p, r < s approximations we require a relationship between these solutions 
to be satisfied at the (ps) order. This relationship does not seem to be manifestly 
satisfied. We therefore conclude that the condition (|4.8|) is a test of the double series 



method itself. If (4.8) is not satisfied at any particular level of approximation the 
expansion procedure does not seem to be consistent with our premise of having an 
isolated source. We now examine the work of Bonnor and Rotenberg (J] on the double 
series method and show that the solutions obtained are consistent in the above sense. 

5 The First and Second Order Solutions. 

We now briefly review the first and second order (in m) solutions obtained by the 
double series method, paying particular attention to the condition (fhq) derived in 
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the previous section. The first and second order solutions were found and interpreted 
by Bonnor and Rotenberg Q] . 

The linear approximation, (Is), is a superposition of terms each involving one of 
the multipolc moments, Q n , about the axis of symmetry. We use the non-gravitational 
forces inherent in the system to ensure that these have a suitable form: 

Q n = ma n h n (u) (5.1) 

for n > 1. The dipole moment can be eliminated by a suitable choice of frame. 

(ii) 

We choose such a frame and then are able to specify that g a fc = 0. The functions 
h n describe the oscillation of the source. Here, for simplicity we assume that only 
/12 is not identically zero so that the source has only a quadrupole moment. This 

(ps) (ps) (ps) (ps) 

ensures that H , K , N , D all have only even order dependence on the Legendre 

(p«) 

polynomials in cos 6 and thay L depends only on the derivatives of even order 
Legendre polynomials in cos 9. As a result the condition mentioned in the previous 
section for the avoidance of singularities in B will always be trivially satisfied in 
this work. The function h,2 is assumed to be constant for u < u\, u > U2 and to 
vary smoothly for u\ < u < 112- In the double series method all source terms are 
inserted at the linear approximation. This is a distinct advantage over the Bondi, 
van den Burg and Metzner formalism, (||), in which the news function that generates 
the solution contains non-linear terms. To ensure that we do not insert higher order 
source terms we specify that any arbitrary functions of integration at higher orders 
are set to zero, except where these are required to avoid a singularity occurring along 
the axis of symmetry. 

Bonnor and Rotenberg then solved the (22) approximation. The non-linear terms 
that arise at the m 2 a 2 level are due to 

(10) (12) 

gab x g ab (5.2) 
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products. The (22) solution thus describes the interaction between the Schwarzschild 
source and the quadrupole wave. The solution contains integrals of the form: 



5,, 



w~ n h 2 (u + 2r- 2w)dw. 



(5.3) 



These represent wave-tails, describing the back-scattering of the quadrupole radiation 
by the source. 

The (24) solution was solved to 0(r -3 ) by Bonnor and Rotenberg and completed 
by Hunter and Rotenberg Q . It is here that the loss of mass of the source due to the 
gravitational radiation is first seen. The occurrence of such integrals as: 



Y := 



... 2 

h 2 du 



(5.4) 



which are zero for u < u\ and non-zero positive constants for u > u 2 indicate a 
permanent change in the metric due to the period of oscillation. It can be shown that 
they represent a source losing mass at a rate equal to that at which energy is radiated 
by gravitational waves. 

The (24) approximation is the first at which the consistency condition ([D]) is not 
trivially satisfied. At the (22) level, qH , oK, oN all vanish. In the (24) approxima- 
tion, Hunter and Rotenberg found that 



oH 
oK 

oN 



fr~*h 2 



5 r n 2 



W~ 2 h 2 +r- 



k hi h\ 



15 



A h 2 



15 



4h\h 2 



2h 2 



2 4r- 7 h\h 2 



15 

3r- 8 ^ 



4 r^h 2 h 2 



4r -3 hi h 2 + r~ 5 



15 
_-6 



1 h 2 h 2 



\h 2 h 2 



-U 2 h 2 



Tlt 2 



fr~ 7 h 2 h 2 + 3r- 8 h 2 2 



(5.5) 



It can easily be checked that these values satisfy the condition. Thus to second 
order the expansion method of the double series method is seen to allow a consistent 
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specification of an isolated source. It should be noted that the non-linear terms ( p7^ ) 
arise only from the product 



At the third and fourth orders it is not this simple - the non-linear terms consist of 
several such products. If the double series method were to be 'clean', each of these 
products could be considered independently. This does not turn out to be the case. 

6 Third Order Solutions. 

The (32) solution is given in another paper |@). Here it is sufficient to note that at 
this order the condition ( |4.8| ) is trivially satisfied. Here we are interested in the (34) 
solution. 

The non-linear terms of order m 3 a 4 are formed by the following products: 



Let us consider each case in turn. Firstly we consider those non-linear terms that 
arise through (10) x (24) products. Using the power series package developed for 
SHEEP |?J we can calculate the corresponding values of qH , qK and qN. We do not 
give these values here since they are not of particular interest in themselves. When 
we test the consistency condition we find that it is not satisfied. The value of X in 
(fh^) is non-zero. We denote this value X\\ 



(12) (12) 
9ab X g a b ■ 



(5.6) 



(10) (24) (12) (22) (10) (12) (12) 

g a b X gab, g a b x g ab , g ab x g ab x g ab . 



(6.1) 



x 1 



^r~ 4 h 2 hi + r- 6 



fh 2 h 2 + f 



h 2 h 2 



+ 



2 ' 



9 -fr~ 8 h 2 h 2 - 



168r -9 ^. 



(6.2) 



36h 2 h 2 - A8h 2 
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Now we consider the (12) x (22) case. Again we find that the consistency condition 
fails. We denote the values of X that we find in this case X^. 



X 9 



2 ~- 4 h 2 h\ 



16 „ 



IS 



hi hi 



— tl 2 tl 2 



-I2r- 7 ti 2 h 2 + ^r~ 8 h 2 h 2 + ^r^hi. 



(6.3) 



Finally we look at the non-linear terms that arise through (10) x (12) x (12) products. 
We find a third, non-zero value for X which we this time denote X 3 : 



Xo 



4r" 4 h 2 h 2 + fr- 5 h 



-24h 2 h 2 + 48ft,. 



hi hi 



344 h h 



5 



504 -9 1 2 

— r n 2 . 



(6.4) 



Thus we find that each individual product in (6.1) fails to satisfy the condition 
( [i~§| ) and hence leads to singularities in the metric. However, if we consider all the 
products together, i.e. the complete (34) case, we find that the condition is satisfied. 
This is seen by the relationship: 



x 1 + x 2 + x 3 = 0. 



(6.5) 



The singular terms that appear due to the non-zero values of X given by ([T^) to 
(^J) exactly cancel each other. Thus the double series method to this order is mixed 
in the sense that the decomposition of the problem into the individual cases given by 
( |6.l| ) leads to a solution that in each case is inconsistent with our assumption that the 
source is isolated. Only the full solution is valid. We have obtained the solution to 
the (34) field equations in order to proceed to the fourth order solution but we do not 
give it here. The fact that it cannot be decomposed into its constituent cases would 
lead to great difficulties in making a physical interpretation of the solution. This is 
certainly the case at the (44) level which we consider next. 
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7 Fourth Order Solutions. 



Here we consider the (44) solution. This would appear to be of significant physical 
interest since it contains the first interaction between the wave-tails, arising from the 
product 



(22) (22) 

g a b x g a b 



(7.1) 



However, when we calculate the non-linear terms due to this product we find that the 



consistency condition (4.8) is not satisfied. In particular 



X 



4„-3 



ho Sk 



h 2 ^4 — ^h 2 S 5 
2 



5 h 2 5*3 — 8/1264 + ^-h 2 S^ 



-| h 2 h 2 + § h 2 S 2 + ±h 2 - W h ^ S 3 + 9/12^4 + f h 2 S 5 



6h 2 h 2 - Ylh 2 S 2 + ^-h 2 S; 



18 u c 312 

-T-n 2 b 2 K - 



h 2 S 3 



5 '-2»J3 
-9 



^■h 2 S4 



3h 2 h 



2 "2 



'-fh 2 h 2 



^h 2 S 2 



b 4r- w hl 



(7.2) 



Thus the (22) x (22) case leads to singularities along the axis appearing in the 



metric at the 



level. 



However, just as at the third order, we should also consider all possible products 
that arise in the non-linear terms: 



(10) (34) (12) (32) (22) (22) 

g a b x g a b, g a b x g a b, g a b x g a b 



(7.3) 



as well as triple and quadruple products. We then find that the condition (4 



again satisfied. The double series method allows a consistent definition of an isolated 
source to fourth order. The singular terms that do appear in the wave-tail - wave-tail 
interaction ( ]T. 2| ) are exactly cancelled in the full (44) solution. This suggests that it is 
not possible to separate the physically interesting part of the (44) solution from that 
part with less apparent physical interest. This would make any physical interpretation 
of the (44) solution very difficult. 
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8 Discussion. 



We have derived a condition for the double series method that must be satisfied 
at each order for the solution to that order to be consistent with the idea of an 
isolated source. Failure to satisfy this condition leads to a singularity along the axis of 
symmetry at that order. We have then verified that this condition is satisfied to order 
m 4 a 4 for a source with only a quadrupole moment. This certainly suggests that the 



particular relationship (4.8) is a feature of the field equations and that the condition 
is automatically satisfied when we are dealing with isolated sources. However it seems 
very difficult to prove this. It would appear likely that any generalisation of the double 
series method to isolated sources without axisymmetry would require an analogous 
condition to be satisfied. 

When we considered the third and fourth order solutions we found that the most 
straight-forward decomposition of the solution into a series of independent physical 
effects proves to be impossible. There is a complicated mixing of the individual 
interactions that ensures that despite the fact that the separate cases each seem 
to generate singular terms in the metric, the full solution remains non-singular and 
consistent with the definition of a source that remains isolated. The cancellation 
of these singular terms that occurs in the full solution at each order, is certainly 
surprising. It is also very important for the validity of the expansion procedure 
inherent in the double series method. 
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A Appendix 

Here we present the vacuum equations used in the double series method and give 
details of solving these equations. To order m p a s the vacuum Einstein equations are 
given by the vanishing of the Ricci tensor components which are (dropping the labels 

2Ru = -Ar~ 1 F 1 - H 
2r~ 2 i? 22 = Bu -2B li + 2r- 1 {B 1 -B i + D 1 - F 1 - G 12 ) + 

r~ 2 (-B 22 - 'iB 2 cot 9 + 2B + 2D + 2F 22 - AF - AG 2 - 
2Gcot<9) - I 

-2R\ = -B n +2B 1A + 2r~ 1 {-B 1 + B 4 + D 1 - Fx - dcote) + 
r~ 2 (-B 22 - 3B 2 cot 6 + 2B + 2D + 2F 2 cot 9 - 4F - 
2G 2 - AG cot 9) - J 
2R 44 = -D n + 2F 14 + 2r- 1 (-D 1 - D 4 + 2F 4 + G 24 + G 4 cot 9) - 
r- 2 (D 22 +D 2 cot8)-K 
2r- 1 R l2 = -G 11 +r- 1 (-B 12 -2B 1 cot9 + F 12 -2G 1 ) + 
2r- 2 (-F 2 + G)-L 
2R U = -D 11 + 2F 14 + r- 1 (-2D 1 + G 12 + G 1 cot6) + 
r~ 2 (-F 22 - F 2 cot 9 + G 2 + G cot 9) -N 
2r- x R 24 = -G 11 + G 14 + r- 1 {-B 24 -2B 4 cot9-D 12 + F 12 +F 24 - 

2G 1 -G 4 )-P (A.l) 

where H, I, J, K, L, N, P are the non-linear parts of the equations (the right hand 



sides of (2.6)) 



17 



The first equation is directly integrable with respect to r to give: 



F = -\ J rHdr + rj(9,u) 



(A.2) 



where rj is a function of integration. We can also integrate the R14 component with 
respect to r to give: 

r(G 2 +G cot 9)=r 2 D 1 + j {r 2 {N - 2F 1A ) + (F 22 + F 2 cot 6)}dr + X (9, u) (A.3) 
where X is a l so a function of integration. We can now use this to eliminate G from 



the i?44 component. This leads to the pseudo-wave equation (3.2): 

UD = D n -2D li + 2r- 1 (D 1 +D i )+r- 2 (D 22 +D 2 cot( 
= -K + 2(F li + 2r- 1 F 4 ) + 



2r" 



J{r 2 (N - 2F U ) + {F 22 + F 2 cot6)}dr + X 



(A.4) 



To hnd G we integrate (A.3) with respect to 



G = r^ 1 J F 2 dr + r -1 cosec0 / sin 9 J ' r 2 (N - 2F lA )dr + r 2 D x + x 
+v(r, u) cosec^ 



dB 



(A-5) 



with v being another function of integration. Finally we integrate the R\ 2 component 
with respect to r and 9 to obtain an expression for B: 



B = cosec 2 9 J sin 2 9 - J {rL + 2r~ 1 (F 2 - G)}dr + F 2 - G - rG x 
+r(r, u) cosec 2 9 + [i(9, u) 



d9 



(A.6) 



where r, [i are two further functions of integration. 
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